CACHE (Critical Assessment of Computational Hit-finding Experiments): A public–private partnership benchmarking initiative to enable the development of computational methods for hit-finding

One aspirational goal of computational chemistry is to predict potent and drug-like binders for any protein, such that only those that bind are synthesized. In this Roadmap, we describe the launch of Critical Assessment of Computational Hit-finding Experiments (CACHE), a public benchmarking project to compare and improve small-molecule hit-finding algorithms through cycles of prediction and experimental testing. Participants will predict small-molecule binders for new and biologically relevant protein targets representing different prediction scenarios. Predicted compounds will be tested rigorously in an experimental hub, and all predicted binders as well as all experimental screening data, including the chemical structures of experimentally tested compounds, will be made publicly available and not subject to any intellectual property restrictions. The ability of a range of computational approaches to find novel binders will be evaluated, compared and openly published. CACHE will launch three new benchmarking exercises every year. The outcomes will be better prediction methods, new small-molecule binders for target proteins of importance for fundamental biology or drug discovery and a major technological step towards achieving the goal of Target 2035, a global initiative to identify pharmacological probes for all human proteins. [Figure not available: see fulltext.

Quantifying Water-Mediated Protein–Ligand Interactions in a Glutamate Receptor: A DFT Study

It is becoming increasingly clear that careful treatment of water molecules in ligand–protein interactions is required in many cases if the correct binding pose is to be identified in molecular docking. Water can form complex bridging networks and can play a critical role in dictating the binding mode of ligands. A particularly striking example of this can be found in the ionotropic glutamate receptors. Despite possessing similar chemical moieties, crystal structures of glutamate and α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) in complex with the ligand-binding core of the GluA2 ionotropic glutamate receptor revealed, contrary to all expectation, two distinct modes of binding. The difference appears to be related to the position of water molecules within the binding pocket. However, it is unclear exactly what governs the preference for water molecules to occupy a particular site in any one binding mode. In this work we use density functional theory (DFT) calculations to investigate the interaction energies and polarization effects of the various components of the binding pocket. Our results show (i) the energetics of a key water molecule are more favorable for the site found in the glutamate-bound mode compared to the alternative site observed in the AMPA-bound mode, (ii) polarization effects are important for glutamate but less so for AMPA, (iii) ligand–system interaction energies alone can predict the correct binding mode for glutamate, but for AMPA alternative modes of binding have similar interaction energies, and (iv) the internal energy is a significant factor for AMPA but not for glutamate. We discuss the results within the broader context of rational drug-design

Sperm performance in conspecific and heterospecific female fluid

Divergent sexual selection within allopatric populations may result in divergent sexual phenotypes, which can act as reproductive barriers between populations upon secondary contact. This hypothesis has been most tested on traits involved in precopulatory sexual selection, with less work focusing on traits that act after copulation and before fertilization (i.e., postcopulatory prezygotic traits), particularly in internally fertilizing vertebrates. However, postcopulatory sexual selection within species can also drive trait divergence, resulting in reduced performance of heterospecific sperm within the female reproductive tract. Such incompatibilities, arising as a by-product of divergent postcopulatory sexual selection in allopatry, can represent reproductive barriers, analogous to species-assortative mating preferences. Here, we tested for postcopulatory prezygotic reproductive barriers between three pairs of taxa with diverged sperm phenotypes and moderate-to-high opportunity for postcopulatory sexual selection (barn swallows Hirundo rustica versus sand martins Riparia riparia, two subspecies of bluethroats, Luscinia svecica svecica versus L. s. namnetum, and great tits Parus major versus blue tits Cyanistes caeruleus). We tested sperm swimming performance in fluid from the outer reproductive tract of females, because the greatest reduction in sperm number in birds occurs as sperm swim across the vagina. Contrary to our expectations, sperm swam equally well in fluid from conspecific and heterospecific females, suggesting that postcopulatory prezygotic barriers do not act between these taxon pairs, at this stage between copulation and fertilization. We therefore suggest that divergence in sperm phenotypes in allopatry is insufficient to cause widespread postcopulatory prezygotic barriers in the form of impaired sperm swimming performance in passerine birds

Allosteric mechanism in AMPA receptors: A FRET-based investigation of conformational changes

α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors are the primary mediators of fast excitatory synaptic transmission in the mammalian CNS. Structures of the extracellular ligand-binding domain suggest that the extent of cleft closure in the ligand-binding domain controls the extent of activation of the receptor. Here we have developed a fluorescence resonance energy transfer-based probe that allows us to study the extent of cleft closure in the isolated ligand-binding domain in solution. These investigations show that the wild-type protein exhibits a graded cleft closure that correlates to the extent of activation, which is in qualitative agreement with the crystal structures. However, the changes in extent of cleft closure between the apo and agonist-bound states are smaller than that observed in the crystal structures. We have also used this method to study the L650T mutant and show that in solution the α-amino-5-methyl-3-hydroxy-4-isoxazole propionate-bound form of this mutant exists primarily in a conformation that is more closed than predicted based on the activity, indicating that the degree of cleft closure alone cannot be used as a measure of extent of activation of the receptor, and there are possibly other mechanisms in addition to cleft closure that mediate the subtleties in extent of activation by a given agonist